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ABSTRACT 



Context. We study the X-ray emission of the Galactic supernova remnant (SNR) CTB 109 (G109. 1-1.0), which is well-known for its 
enigmatic half-shell morphology both in radio and in X-rays and is associated with the anomalous X-ray pulsar (AXP) IE 2259+586. 
Aims. We want to understand the origin of the X-ray bright feature inside the SNR called the Lobe and the details of the interaction 
of the SNR shock wave with the ambient interstellar medium (ISM). 

Methods. The Lobe and the northeastern part of the SNR were observed with Chandra AOS -I. We analysed the spectrum of the X-ray 
emission by dividing the entire observed emission into small regions. The X-ray emission is best reproduced with one-component or 
two-component non-equilibrium ionisation models depending on the position. In the two-component model, one emission component 
represents the shocked ISM and the other the shocked ejecta. 

Results. We detect enhanced element abundances, in particular for Si and Fe, in and around the Lobe. There is one particular region 
next to the Lobe with a high Si abundance of 3.3 (2.6 - 4.0) times the solar value. This is the first, unequivocal detection of ejecta in 
CTB 109. 

Conclusions. The new Chandra data confirm that the Lobe was created by the interaction of the SNR shock and the supernova 
ejecta with a dense and inhomogeneous medium in the environment of SNR CTB 109. The newly calculated age of the SNR is 
1.4xl0 4 yr. 

Key words. Shock waves - ISM: supernova remnants - X-rays: ISM - X-rays: individual: SNR CTB 109 
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1. Introduction 

The Galactic supernova remnant CTB 109 (G109. 1-1.0) is 
the host of the anomalous X-ray pulsar (AXP) IE 2259+586 
(Fahlman & Gregorylll98ll) and represents one of the most ex- 
otic and interesting objects in the X-ray sky. The study of an 
SNR associated with an AXP provides valuable information on 
the environment in which it formed and an independent estimate 
of the age of the objects. Given that there are on ly three firm as- 
sociations of SNRs with AXPs (lMereghettil2008l) . detailed stud- 
ies of each association will produce progress in our understand- 
ing not only of SNRs but also of the AXPs. 

The SNR CTB 109 has a spectacular semi-circular morphol- 
ogy in both the X-ray an d the radio. At a distance of 3.2+0.2 kpc 
(iKothes & FosteJ 120121) , it is located next to a giant molecu- 
lar cloud (GMC) complex and is one of the most striking ex- 
amples of an interaction of an SNR with a molecular cloud . 
Since there is neith er X-ray (see Fig. 1 in ISasaki et aD 120041) 
nor radio emission (iHughes et al.lll98lT) where the western part 
of the shell would be expected, the semi-circular morphology 
of CTB 109 implies that the shock has been greatly impeded 
or even stopped entirely by the GMC complex in the west. 
A linear feature in CO ('CO arm') extends from the GMC 
complex t o the local X-ray min imum in the northern half of 
the SNR dTatematsu et al.l fl987l) . Therefore, at least a part of 
the GMC complex extends in front of the remnant. CTB 109 
has an X-ray bright interior region known as the 'Lobe'. The 
Lobe is brighter than any part of the shell. Although it has 



been suggested tha t this feature is a jet ass ociated with the 
AXP IE 2259+586 (Gregory & Fahlm anll 1983b high -resolution 
images from RO SAT HRI (iHurford & Fesenl fl995h. Chandr a 
dPatel et al.ll200ll) . as well as XMM-Newton (ISasaki et al.ll2004 . 
show no morphological connection with the pulsar. Furthermore, 
the X-ray spectrum fr om the Lobe obtain ed with XMM-Newton 
is completely thermal (ISasaki et al .1120041) . 

We believed that the bright X-ray emission is the result of the 
interaction between the SNR shock and a molecular cloud com- 
plex and therefore p erformed an observati on with the Chandra 
X-ray Observatory (iWeisskopf et al.l|2002l) using the Ad vanced 
CCD Imaging Spectrometer (ACIS. iGarmire et al.ll2003l) in or- 
der to study this interaction region in more detail. The first anal- 
ysis of the high resolution data obtained with ACIS-I in combi- 
nation with new high resolution CO data from the Five College 
Radio Observatory has revealed regions with signs of interac tion 
between the SNR shock and CO clouds ( Sasa ki et al.l 200$. In 
this paper we present the results of the spectral analysis of the 
entire northeastern part of the SNR, which was observed with 
ACIS-I. 

2. Chandra data 

We observed the northeast part of the SNR in an 80 ksec point- 
ing of Chandra with ACIS-I as the prime instrument (ObsID 
4626). The Lobe and the northeastern part of the remant shell 
were completely covered by the ACIS-I array (Fig.[T}. The data 
were analysed with CIAO 4.4. and CALDB 4.4.8. 
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Fig. 1. Exposure-corrected Chandra ACIS-I intensity map (0.35 
-8.0 keV). The big red circle indicates the position of the Lobe, 
the yellow ellipse that of the CO arm. 



2.1. Images 

Figure [T] shows the exposure-corrected intensity map of the 
ACIS-I data (0.35 - 8.0 keV). The events were binned with a 
bin size of 4 pixels and the image has been smoothed with a 
gaussian kernel of 3 and 4 pixels in Figs.[T]and[2](see below), 
respectively. 

Compared to the XMM-Newton image dSasaki et al.ll2004l) . 
the Chandra image reveals more point sources. Also, structures 
in the dark region absorbed by the CO arm, as well as filamen- 
tary structures in the Lobe, are resolved. The two-colour image 
(Fig.0), in which the soft band (0.35 - 0.8 keV) is presented in 
red and the hard band (1.0 - 8.0 keV) in green, shows colour 
variations in the Lobe. There is a bright yellow structure in the 
central part of the Lobe, while the outer diffuser parts are green 
(harder) in the east and red (softer) in the west. The two-colour 
image shows spectral variations that are either intrinsic or due 
to absorption or perhaps a combination of the two. The diffuse 
emission in the CO arm appears to have harder spectra consistent 
with higher absorption. A detailed spectral analysis would possi- 
bly distinguish between intrinsic spectral variations and variable 
absorption. 

For further analysis, we divided the diffuse emission of the 
SNR into small regions with similar surface brightness and X- 
ray colour, and extract spectra for each region. To study the spec- 
tra of only the diffuse emission, we first performed source detec- 
tion on the whole data (wavdetect) and excluded all detected 
point and point-like sources in each region. The regions used are 
shown in Figure [3] The number of counts in the regions varies 
from ~4 000 - 20000 cts. The spectra are binned with a min- 
imum of 20 counts per bin so that Gaussian statistics may be 
assumed in the fitting. 
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Fig. 2. Two-colour image with red (0.35 - 0.8 keV) and green 
(1.0-8.0keV). 
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Fig. 3. Extracted regions overlaid on the intensity map. Point 
sources are removed. 

2.2. Spectral analysis 

We analysed the spectra using the X-ray spectral fitting package 
XSPEC Ver. 12.7.1. The analysis of XMM-Newton data showed 
that the emission from C TB 109 is thermal with no indication 
of nonthermal emission ( Sasaki et al. 2004). Furthermore, the 
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Fig. 4. Contour diagram of the parameters kT and r for the fit 
with one VNEI model of region 49. The contours correspond to 
confidence levels of 60%, 90%, and 99%. 



thermal emission suggests that the plasma is not in collisional 
ionisation equilibrium (CIE). Therefore, we fitted the extracted 
ACIS-I spectra with non-equilibrium ionisation (NEI) models. 
The foregr ound absorbing col umn density Nu was modeled us- 
ing TBABS (Wil ms et ^Dl20QQh. We used ab undances relative to 
solar values reported by lWilms et al.l fcOOO). 

2.2.1. One-component NEI model 

First, we fitted all spectra assuming one VNEI component, 
similar to what w as used to analyse the XMM-Newton data 
(iSasaki et al.ll2004l) . Initially, the abundances were all fixed to 
solar values. In regions 49 and 57, which correspond to the 
outer shock, we get relatively good fits with Nu = 9.6(8.8 - 
10.5flx 10 21 cm" 2 , 1.2(1.1 - 1.3) x 10 22 cm" 2 , kT = 0.25 (0.22 
- 0.27) keV, 0.25 (0.23 - 0.26) keV, and n e t > 7.2 x 10 11 s cm" 3 , 
> 2.3 x 10 10 s cm" 3 , with reduced^ 2 = 1.4 for 107 and 99 
degrees of freedom for regions 49 and 57, respectively (see 
Fig. |4]) . The best fit values for the ionisation timescale are n Q t = 
4.9 x 10 13 s cm" 3 and 2.1 x 10 13 s cm" 3 indicating CIE in these 
regions. 

We also freed the abundances, but the one-component de- 
scription of the data is not satisfactory in about 25% of the re- 
gions with reduced^ 2 > 1.5, ranging up to reduced^ 2 = 2.0 (see 
Fig. [51 upper left diagram). 

2.2.2. Two-component NEI model 

Therefore, in the next step we assumed two thermal VNEI com- 
ponents: 1) to model the emission of the shocked ISM and 2) to 
verify the existence of emission from shocked ejecta. This two- 
component model improves the fits in a number of regions (see 
below), especially in region 39, which had a reduced x 2 =2.0 
(103 degrees of freedom) for the one-component VNEI model 
and can be fitted with a reduced^ 2 = 1.1 (100 degrees of free- 
dom) with a two-component VNEI model (Fig.[5j left panels). 

For the ISM component, we use kT\ =0.25 keV and n Q t\ = 
1 x 10 12 s cm" 3 as starting values, representing the best fit pa- 
rameters of the single-component VNEI model for regions 49 
and 57. All abundances of this component were fixed to solar 

1 All errors in this paper are 90% confidence errors. 



values. The second component, which was introduced to model 
the ejecta emission, was fitted with variable abundances for ele- 
ments showing strong emission lines. The fits required a higher 
temperature of kT^ ~ 0.6 keV. The ejecta component with the 
higher temperature dominates the spectrum for energies higher 
than ~1 keV in most of the spectra (see Fig. [3 middle panels). 
The abundances of Mg, Si, S, for which emission lines are visi- 
ble in the spectra, as well as Fe, were freed for the second (ejecta) 
model component and fitted. In some regions, higher abundances 
are found for Si and Fe in particular (see Fig. [5] right panels). 

2.2.3. F-test 

To verify how much the fits improve by including a second VNEI 
component, we performed an F-test for the one-component VNEI 
and two-component VNEI models for all regions. In 27 out of 
58 regions, the two models differ with probabilities higher than 
99.9%. In order to visualize the variations of the spectral prop- 
erties, we created a map of the spectral parameters by filling 
each region with the corresponding parameter value. In the top 
row of Fig. [6] we show the distribution of the temperature pa- 
rameter kT of the one-component VNEI model for regions in 
which the F-test indicated no necessity of an additional com- 
ponent (left) and the distributions of kT\ (ISM) and kT^ (ejecta) 
of the two-component VNEI model (middle and right, respec- 
tively). Regions in which the F-test confirmed an improvement 
of the fits after including the second VNEI component are marked 
with white crosses. The middle and bottom rows in Fig. show 
the same for the ionisation timescales r (left) and T2 (right) and 
the foreground column densities Nu for the one-component (left) 
and two-component (right) fits. 



3. Discussion 

The analysis of the Chandra ACIS spectra of the northeastern 
part of the SNR CTB 109, divided into small regions, has shown 
that the spectra can be well fitted by a one-component VNEI 
model in many regions. However, in 50% of the regions a two- 
component model consisting of two thermal VNEI models with 
two different temperatures improves the fit significantly. In this 
model, the first component is used to describe the emission from 
the shocked ISM, while the second component reproduces the 
emission from the ejecta. Therefore, we believe that there are at 
least two emission components in all regions. In cases in which 
the one-component model yields a sufficiently good fit, the tem- 
perature and the ionisation timescale most likely tend toward 
values of the component of the multi-temperature spectrum that 
dominates the emission. We base our discussion on the results 
obtained with the two-component spectral fits for all regions. 
The diagrams in Fig. [7] show the relations between some param- 
eters derived from the fits (see following subsections for further 
discussion). 

The temperature of the ISM component (Fig. [6) upper mid- 
dle) is comparable in all regions and is low (~ 0.1 - 0.3 keV), 
while the temperature of the ejecta component is higher (~ 
0.4-0.9 keV, Fig.[6j upper right). The ionisation timescale of the 
ejecta component is T2 = 10 11 " 12 scm" 3 (Fig.0 middle panels) 
except for regions at and around the CO arm, in which T2 tends 
to be higher. The higher values for T2 in the CO arm regions in- 
dicate a higher density as we can assume that these regions have 
been shocked at a similar time to the rest of the remnant. 

The Nu images (Fig.[6l lower panels) show that the fore- 
ground absorption is high near the CO arm (Nu> 1.1 X 
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Chandra ACIS-I Region 39 (VNEI) Chandra ACIS-I Region 57 Chandra ACIS-I Region 17 
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Fig. 5. Left: Chandra ACIS-I spectrum of an interior region 39 with the best fit single-component VNEI model (upper diagram, 
reduced^ 2 = 2.0 for 103 degrees of freedom) and two-component model (lower diagram, reduced^ 2 = 1.1 for 100 degrees of 
freedom). Middle: Spectrum of the outermost region 57 in the east (upper) and an interior region 54 (lower) with the best fit 2 VNEI 
model. The relative flux of the ejecta component (dotted) with respect to the ISM component (dashed) in region 57 is lower than in 
most other regions. Right: Spectrum of region 17 with enhanced Si xm emission (upper diagram). The Mg xi, Si xm, and S xv triplets 
are marked. For comparison, the spectrum of region 08, which is the brightest region in the Lobe next to region 17, is shown (lower 
diagram). The best fit 2 VNEI model is additionally plotted in the diagrams. 



10 22 cm" 2 ). Around and, in particular, west of the Lobe, the fore- 
ground absorption is lower (Nu< 1 .Ox 10 22 cm" 2 ), while it seems 
to be higher southeast of the Lobe. 



value distribution map in Fig. [8] in which the inner (western) re- 
gions tend to have higher fz values than the outer regions. 



3.1. Filling factors and ejecta mass 

CTB 109 is an evolved SNR in the Sedov phase. In such a 
remnant, which also shows interactions with denser molecular 
clouds, we can assume that the reverse shock has already propa- 
gated to the centre and thus heated most or all of the ejecta. The 
spectral model provides us with the normalisation parameter 
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with f\ and fz being the filling factors of the ISM and ejecta 
components, respectively (f\ + fz = 1). For the di stance to the 
SNR, we use D = 3.0 ± 0.5 kpc as esti mated b y iKothes et al.1 
(120021) and confirmed by Kothes & Foster (I2012L 3.2 ± 0.2 kpc). 
To calculate the filling factor for the ejecta component fz we 
assume pressure equilibrium between the two thermal compo- 
nents. Applying momentum conservation for the shock propa- 
gating through these components, we get 



[norm\ T\ 
normz T\ 



+ 1 



-1 



(2) 



fsee lBocchino et al.lll999h . 

The resulting parameter map for the filling factor of the 
ejecta component of the two-component VNEI fits is shown in 
Fig. [8] for all regions. We indeed see a trend in the parameter 



We can estimate the total ejected mass using the filling fac- 
tor fz and the density nz derived for the ejecta component from 
Eq.Q] We assume that the SNR is a half sphere with a radius of 
18 f5 + lfO as d etermined from the XMM-Newton EPIC images 
bv lSasakietaD (120041) . For D = 3.0 ± 0.5 kpc, R = 16 ± 3 pc 
= 5.0 ± 0.8 x 10 19 cm. Since we have performed spectral analy- 
sis using a model with two thermal components only in regions 
that have been observed with Chandra ACIS-I, we use the me- 
dian values of fz and nz from the two-component fits to esti- 
mate the total mass of the ejecta. With /2,median = 0-44 and 
^2,median = 0-62 cm" 2 , the total ejecta mass under the assump- 
tion of a half sphere as seen, e.g., in the XMM-Newton mosaic 
image, is M e j ecta = 60 Mq. The density nz in the regions ob- 
served with Chandra varies from 0.23 to 0.81 cm" 2 , while the 
filling factor fz has values between 0.29 and 0.91. Therefore, 
these estimates have uncertainties of -50%, whereas the volume 
may also have an uncertainty of the same order since the X-ray 
emitting hot gas might not fill the entire volume applied for the 
estimate. In addition the ejecta might be clumped and thus not 
uniformly distributed in the assumed volume. If we suppose that 
the ejecta only fills 30% of the volume, the ejecta mass will be as 
low as M e j ecta = 20Mq. Therefore, one should keep in mind that 
Ejecta is a very crude estimate. However, even the lower limit of 
Mejecta = 20Mq would rule out a Type la supernova as the origin 
of CTB 109. 
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Fig. 7. Distribution of the parameter values for the fit with two 
VNEI components. The labels of the data points are region num- 
bers. The lines in the diagrams for surface brightness norm\^/A 
against temperature kT\^ (upper and middle panels) show the 
isobars for p/k = 10 6 cm -3 K (blue) and 10 7 cm -3 K (green) for 
different values of the filling factors in the remnant and the size 
of the emitting volume along the line of sight. The line for each 
filling factor with the lowest norm\^/A are marked with a thicker 
line and labelled with the f\ t 2 value. The lines in the diagram 
showing the surface brightness as normal A against the ionisa- 
tion timescale T2 (bottom panel) are isochrones for t - 10000, 
20000, 50000, and 100000 yr. A range for all isochrones is 
shown for each t covering the range of variation in the depth 
of the emitting volume. 




Fig. 8. Distribution of the filling factor for the ejecta component 
fz derived from the fits with the two-component VNEI model. 



3.2. Pressure 

In the norm\ r 2 1 'A-kl \ t 2 diagrams in Fig. [7] (upper two dia- 
grams) we have plotted lines of constant pressure with p/k = 
10 6 cm" 3 K (blue) and 10 7 cm" 3 K (green) for different values 
of filling factors and depths of the emitting volume in the SNR. 
These lines were calculated from the best fit parameters assum- 
ing pressure equilibrium between the ISM and ejecta compo- 
nents. All the data points are consistent with a thermal pressure 
of p/k « 10 7 cm" 3 K, corresponding to 1.4 xlO" 13 dyne cm" 2 . 
The scatter is mainly due to the different filling factors and the 
uncertainties in the size of the emitting volume along the line of 
sight. 



3.3. Age estimate 

In the T2-norrri2/A plot in Fig. [7] (lower diagram), we plot the 
ionisation timescale T2 of the ejecta component versus the nor- 
malisation norrri2 per area A of each region. 

To study the dependence of the X-ray surface brightness, 
measured as the normalisation parameter per area norrri2/A, and 
the ionisation timescale T2, we can replace n Q by r/t in Eq.[T] 
We then get norrri2/A as a function of r for a given t. The value 
of norm2 depends on the emitting volume, thus on the depth of 
the volume along the line of sight for each region. Assuming that 
the AXP is located at the centre of the circle corresponding to the 
SNR shell, we calculated the minimum and maximum distances 
from the AXP for each region. The calculated isochrones for dif- 
ferent t values are also shown in the T2-norm2/A plot. The two 
lines for each time t show the possible ranges for the depth of 
the remnant. This diagram indicates that the age of the shocked 
plasma, derived from the ejecta component (2) is ~ 2 x 10 4 yr. 

The fits have shown that the temperature of the ISM compo- 
nent (1) is between 0.1 keV and 0.3 keV. For the outer regions 
49 and 57, the single component VNEI fit yields kT = 0.25 keV 
(Sect. l2.2.T1) . The high and not well-constrained values of the 
ionisation timescale r in the single-component fits, as well as 
the T\ in the two-component fits, suggest that the plasma of the 
shocked ISM is close to CIE. If we assume that this temperature 
is more or less representative of the shocked plasma, the blast 
wave velocity is 



3 in 
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with a mean mass per free particle for a fully ionised plasma of 
in - 0.61 m p . For a temperature of kT\ = 0.25 + 0.03 keV, we 
et v = 460 ± 30 km s" 1 . The radius is R = 5.0 ± 0.8 x 10 19 - 



cm 



Sasaki etal]|2004 ) . Using the Sedov si milarity solution (iSedovl 
aylorl fl 



1959; Taylor 1950; von Neumann 1947), the age of the remnant 
can be estimated as 



t = 



2R 
5v* 



(4) 



We thus obtain t = (14 ± 2) x 10 3 yr from the ISM component, 
which is a lit tle higher than the value obtained from the XMM- 
Newton data (Sas aki et~aI1l2 004) for which we only assumed one 
spectral component. Interestingly, this new value for the age of 
the SNR is in agreement with the age estimate obtained from the 
fits of the ejecta component. 

3.4. Ejecta 

As can be seen in the Si-abundance and Fe-abundance distribu- 
tion images in Fig.|9l the abundances measured for the ejecta 
emission are all comparable to or lower than solar values, ex- 
cept for the Lobe. While the Ne, Mg, and S abundances seem to 
agree with those of the other elements, Si and Fe abundances are 
higher, especially in and around the Lobe. The enhanced abun- 
dances in the Lobe suggest that its emission has a contribution 
from an ejecta clump or a conglomeration of ejecta clumps and 
shocked clouds. There is one particular region (region 17) with 
significantly enhanced Si abundance within the Lobe (3.3 [2.6 - 
4.0], see Fig.0), c). The spectrum of region 17 is shown together 
with that of an adjacent region with lower Si abundance in Fig. [5] 
(right panels). 

Another possible process that can cause enhanced emission 
for particular elements in thermal plasma is charge exchange. 
However, to produce charge exchange emission there must be a 
phase of cold neutral gas next to or inside the hot X-ray emitting 
plasma. CO data show no emission at the position of the Lobe 
or west of it, where the Si a bundance enhancement is observed 
in X-rays dSasaki et al.ll2006h . For a mature SNR like CTB 109 
we can assume that the reverse shock has propagated all the way 
through the SNR and, therefore, there should be no unionised 
gas in the interior of CTB 109. Furtheremore, if the enhanced 
emission of Si is due to charge exchange, one would expect even 
higher indication of charge exchange for O or Ne for solar or 
cosmic abundances. Since no emission enhancement is observed 
for these lower Z elements, we can rule out that the enhanced 
emission of particular elements is caused by charge exchange. 

The detection of ejecta emission is also very important for 
the study of the AXP IE 2259+586, as it will give information on 
its progenitor. Together with soft gamma-ray repeaters (SGRs) 
and the sub-class of rotation-powered high-magnetic field pul- 
sars, the AXPs are believed to form a class of neutron stars 
with extremely high magnetic fields. So far, there are only a 
few confirmed associations between SNRs and high-magnetic 
field neutron starfi SNR G292.2-0.5 and the pulsar J1119- 
6127, SNR Kes 75 and the pulsar J1846-0258, SNR Kes 73 
and AXP IE 1841-045, SNR CTB 109 and AXP IE 2259+586, 
SNR G327.2-0.1 and AXP IE 1547.0-5408, SNR G337.0-0.1 
and SGR 1627-41, SNR G042.8+00.6 and SGR 1900 +14, and 
SNR N 49 and SGR 0526-66 in the LMC. Recently. IPark et al.1 
(120121) detected ejecta emission in a deep observation of the 
SNR N49 with Chandra (120 ks exposure). However, the el- 
ement abundance ratio derived from the emission seems to be 
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more consistent with a Type la supernova, which would rule out 
the association be tween the SNR and the SGR, and requires fur- 
ther investigation. lLopez et al.l (1201 ll) measured enhanced abun- 
dances of Mg, Si, and S in SNR Kes 73 associated with the AXP 
IE 1841-045, whereas ISafi-Harb & Kurnarl (12012b reported on 
possible abundance enhancement in the X-ray spectra of SNR 
G292.2-0.5 and SNR Kes 73. 

Ejecta emission is also detected in other mid dle-aged rem- 
nants, e.g. in the Galactic SNR G349.7+0.2 (lLazendic et al.l 
l2005h . which is interacting with a molecular cloud, or SNRs 
in the Large Magella nic Cloud (L MC) 0548-70.4 0534-69.9 
dHendrick et al.ll2003l) . and N49B (IPark et al.ll2003l) . Clumping 
of ejecta seems to be common in core-collapse SNRs. Ejecta 
clumps have also been found in evo lved SNRs like the Vela, 
both in the interior (iMiceli et a l. 2008) and outsid e the main shell 
(lAschenbach et al.lfl995l) . [Kifonidis et al.l 02003) have shown in 
two-dimensional simulations that macroscopic mixing occurs in 
core-collapse supernova explosions forming clumps of metals. 
Rayleigh-Taylor instabilities grow in the layers of metals during 
the expansion, and these layers fragmentate into ejecta clumps. 

We created narrow band images for the elements Mg and Si, 
as prominent emission lines of these elements are visible in the 
spectra of some regions (Sect. l2.2t . The continuum emission has 
been subtracted from the line emission. Figure [l0]shows a three- 
colour image with red for th e soft band XMM-Newton image 
(0.3 - 0.9 keV) presented bv lSasaki etaD (12004 . green for Mg 
emission ([1 .25 - 1 .45 keV band] - surrounding continuum), and 
blue for Si emission ([1 .7 - 2. 1 keV band] - continuum). There is 
one small region bright in Si emission, covered by the extraction 
region 17 with enhanced Si line emission (see Sect. 12.21) . The 
Mg emission is distributed broadly; soft emission is bright in 
the Lobe and the small region in the northeast shell (region 41), 
which has a relatively low foreground absorption of Nu= 0.80 
(0. 70- 0.94) xnF cm' 2 . 

Kap lan et~aD d2QQ9h measured the proper motions of the 
AXP IE 2259+586 using Chandra ACIS-S observations in 2000 
and 2006 pointed at the AXP. The suggested initial position at 
which the supernova occured and the AXP was born is north 
of the AXP and might even have been in the area coinciding 
with the CO arm, depending on the real age of the SNR. On the 
other hand, new near-infrared observations of the counterpart of 
the AXP with the Keck II telescope using laser guide star adap- 
tive optics have been performed to measure the proper motion 
of the AXP. These data indicate that the initial position of the 
AXP was closer to the position of X-ray Lobe, approximately 
half way between today's position of the AXP and the Si clump 
seen in blue in Fig.fTOl (Tendulkar et al., in prep.). Both results 
show that the position of the initial supernova explosion was lo- 
cated west to northwest of the Lobe and suggest that the interac- 
tion of the blast wave with a dense molecular cloud produced the 
X-ray Lobe and, at the same time, caused the formation of the 
reverse shock, which then ran into ejecta enhanced in Si or Fe. 
However, one has to keep in mind that the alignment of the Si 
enhanced region and the Lobe might well be a projection effect. 



4. Conclusions 

Through the analysis of the deep high-resolution data of the 
Lobe and the northeastern part of the SNR CTB 109 taken with 
Chandra ACIS-I, we confirm ejecta emission inside the remnant. 
There is an enhancement in element abundances especially in the 
X-ray bright Lobe. We have thus unambiguously detected ejecta 
emission from an SNR associated with an AXP. The ISM and 
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Fig. 10. Three colour image consisting of red image for sof t 
band XMM-Newton image (0.3 - 0.9 keV, ISasaki et alJl2004 . 
green image for Mg (1.25 - 1.45 keV) - continuum, and blue 
image for Si (1.7 - 2.1 keV) - continuum. The Mg and Si im- 
ages were created from the Chandra data. The green dashed box 
indicates the field of view of the Chandra ACIS-I. The yellow 
cross indicates the present position of the AXP IE 2259+586. 

ejecta emission indicate an age of -10 - 20 kyr for the shocked 
X-ray emitting plasma. 

As one can se e in the XMM-Newton intensity image (Fig. 1, 
ISasaki et al.ll2004 the SNR shell is brighter in the north and the 
south, while it is darker with diffuse structures in the outer parts 
of the same sector in which the Lobe is located (east-northeast). 
The analysis of Chandra data has shown that the foreground ab- 
sorption Nu is lower in and around the Lobe (Fig|6l lower pan- 
els), while the Fe L-shell emission is enhanced in the Lobe, with 
one smaller region with significantly enhanced Si abundance. 
The bright feature seen as the Lobe can be explained by the 
following scenario: the blast wave hit a dense cloud, evaporat- 
ing and maybe encompassing it. Due to the interacton with the 
cloud, a reverse shock was formed and propagated towards the 
centre heating the ejecta. The evaporated cloud is now visible as 
the Lobe with s tronger X-ray emission. 

ICastro etall (120121) report the detection of a GeV source with 
Fermi at the position of CTB 109. The Fermi source is located 
southwest of the Lobe where the SNR appears dark, and thus 
most likely highly absorbed in X-rays. Assuming that the y- 
rays are produced in 7r°-decays caused by the interaction of the 
SNR shock with a dense interstellar cloud now visible as the 
Lobe. lCastro et al.l (120121) derive a density of ~ 120 cm -3 for the 
cloud. They point out the discrepancy between this value and 
the d ensity of 0.9 cm" 3 , which we calculated from the X-ray 
data (Sa saki et~aD2006l) . However, as the density of the Lobe of 
0.9 cm" 3 refers to the hot shocked X-ray emitting gas, which is 
likely to have a lower density than t he original cold cloud, the 
new result from ICastro et aD (120121) is consistent with, rather 
than incompatible, with our results. As we have pointed out 



(ISasaki et al.ll2006l) and confirmed in the new spectral analysis, 
the foreground absorption south of the Lobe (regions 38 and 5 1 
in this work) is high with Af H = 1.3(1.1 - 1.5) x 10 22 cm" 2 and 
comparable to that in the CO arm. Therefor e, our new Chandra 
study as well as the new Fermi results by ICastro et aD (1201 2b 
confirm that there was and still is interaction between the shock 
wave of CTB 109 with the ambient dense material. The remain- 
der of this interac tion is still visible as CO clouds presented in 
ISasaki et al.l (12006b . 
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VNEI kT (0.10- 0.90 keV) B^HH VNEI+VNEI kT1 (0.10 - 0.90 keV) WM VNEI+VNEI kT2 (0.10 - 0.90 keV) 
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Fig. 6. Top: Best fit temperature parameters for the one-component VNEI model (kT) in regions in which the F-test indicated that 
the one-component model fits the spectrum sufficiently well (left) and temperature parameters kT\ and kT2 for the two-component 
fits for all regions (middle and right, respectively). Middle: Best fit values for the ionisation timescale r for the one-component VNEI 
model (left) and the ionisation timescale T2 for the dominating ejecta component for the two-component fits (right). Bottom: Best 
fit values for the foreground Nu for the one-component VNEI model (left) and for the two-component fits (right). The crosses mark 
the regions in which the F-test indicated that the two- VNEI component model yields a better fit (see Sect. l2.2.3t . 
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(b) VNEI+VNEI Si (0.0 - 5.0 x solar) 




(C) VNEI+VNEI Si significance (1 -3) 




(d) VNEI+VNEI Fe (0.0 - 5.0 x solar) 




(e) VNEI+VNEI Fe significance (1 - 3 ) 




Fig. 9. Best fit abundances assuming a thermal model consisting of two VNEI components. The shown parameters are: abundances 
for elements other than Si, S, and Fe fitted for the ejecta emission (a), Si abundance fitted for the ejecta emission (b) and its 
significance calculated as (abundance of the Si - abundance of the other elements)/error of the abundance of Si (c), Fe abundance 
and its significance (d and e, respectively). 
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